Dispersive measurements of superconducting qubit coherence with a fast, latching 

readout 
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The "quantronium" is a superconducting qubit consisting of a split Cooper pair box in which a 
large tunnel junction is inserted. This circuit has a special bias point where the Larmor frequency 
is, to first order, insensitive to fluctuations in the bias parameters - the charge of the box island 
and the phase of the large junction. At this optimal working point, the state of the qubit can be 
determined by dispersive measurements that probe the second derivative of the state energy with 
respect to these bias parameters. We use the quantronium phase degree of freedom to perform a 
non-linear, dispersive measurement of its inductive response using bifurcation amplification. This 
novel readout projects the state of the qubit in a few nanoseconds, and its latching property allows 
us to record the resulting information in a few hundred nanoseconds. We have measured, using this 
technique, Rabi oscillations and Ramsey fringes with an improved signal to noise ratio and contrast. 
The speed of this new readout scheme also opens the door for a new class of experiments which 
would characterize the relaxation processes associated with the measurement protocol. 



INTRODUCTION 



Superconducting tunnel junction circuits were first 
proposed for quantum information processing several 
years ago, and at present, are the most advanced solid 
state qubits with the longest measured coherence times 
B II H II II II II • Yet, the physical origin of the noise 
sources limiting coherence are still debated, even though 
the theoretical formalism for treating the effects of noise 
in general is well developed |3, 0, E| • It has been con- 
jectured that impurities or defects found on-chip could 
act as such noise sources lllj. These parasitic elements 
may exist in the junction tunnel barriers, the metallic 
electrodes, the circuit substrate, or in some combination 
thereof. In addition, the shadow-mask evaporation tech- 
nique used to fabricate many superconducting qubits typ- 
ically generates extra electrodynamic resonators in close 
proximity to the qubit junctions |12| |. These resonators 
can have a characteristic frequency comparable to the 
qubit Larmor frequency, and are thus suspected to de- 
cohere the qubit. The precise manner in which a qubit 
interacts with uncontrolled degrees of freedom in its envi- 
ronment depends on the topology of the tunnel junction 
circuit and how information is written to and read from 
the qubit. Circuits which have a high degree of symme- 
try can be significantly decoupled from a noisy environ- 
ment when biased at special operating points. The 
choice of readout scheme is also highly si gnif icant. Dis- 
persive measurements of the qubit state |l4l ITU probe 
the reactive part of the response of the circuit, and are 
thus attractive since they minimally excite the spurious 
degrees of freedom described above. 

We report the first coherence measurements of a su- 
perconducting qubit with a non- linear dispersive readout. 
Our approach involves coupling the "quantronium" qubit 
Pi to a novel superconducting amplifier - the Josephson 
bifurcation amplifier (JBA) 16] . The JBA is based on a 
non-linear electrodynamic resonator with two metastable 
oscillation states [T3. In order to perform a readout, the 
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FIG. 1: Schematic of the measurement setup. The quantron- 
ium qubit is a split Cooper pair box with two small Joseph- 
son junctions in which a large junction is inserted for readout. 
This latter junction is shunted by two capacitors in series and 
forms the non-linear resonator of the JBA readout. The qubit 
state is manipulated by sending pulses to the charge port, 
while readout operation is performed by sending a pulse to 
the phase port and analyzing the phase of reflected signal, 
which carries information on the qubit state. 



resonator is RF-energized to a level where its oscillation 
state now acts as a sensitive pointer of the qubit state. 
This technique does not generate any dissipation on chip 
since the resonator is only damped by circuitry outside 
the chip, i.e. a 5017 transmission line with a matched 
circulator and amplifier, and enables a high-fidelity qubit 
readout with a MHz repetition rate. We have measured 
Rabi oscillations and Ramsey fringes with sufficient speed 
that real time filtering to correct for drifts in the charge 
and flux bias becomes possible. Also, several successive 
readouts may be performed within the energy relaxation 
time of the qubit (Ti). This gives valuable information 
on the readout-induced interaction between the qubit and 
its environment, and accounts for the observed contrast. 
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II. THE HAMILTONIAN OF THE 
QUANTRONIUM QUBIT WITH JBA READOUT 

The principle of our experiment is schematically de- 
picted in Fig. 1 and is based, as discussed above, on the 
quantronium qubit, a three junction circuit which is anal- 
ogous to a one-dimensional atom. The set of three junc- 
tions consists of two small junctions, which we assume to 
be identical and have a Josephson energy comparable to 
the charging energy of the island between them, and a 
large junction, whose Josephson energy is approximately 



100 times larger than that of each small junction. The 
gauge-invariant phase difference of the island with re- 
spect to the mid-point of the capacitance shunting the 
large junction is analogous to the position of the elec- 
tron relative to the nucleus of the atom, while the gauge- 
invariant phase difference 6 across the large junction is 
the absolute position of the nucleus. Neglecting the dissi- 
pation induced in the transmission lines, the total hamil- 
tonian of the split Cooper pair box with a JBA resonator 
is H{t) ^ Hbox{t) + Hres{t) with 
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Here, N and Q/2e are the momenta conjugate to the gen- 
eralized positions 9 and S, respectively. The constants 
Ec, Ej, Ej, C and Cg are the single electron charg- 
ing energy of the island between the small junctions, the 
sum of the Josephson energy of the two small junctions, 
the large junction Josephson energy, the total capaci- 
tance shunting the large junction, and the gate capaci- 
tance, respectively. Here lpq = h/2e is the reduced flux 
quantum. The control parameters U{t) = Urf {t) cosilt 
and I{t) — Irf (t) cos Lut are analogous to electromagnetic 
probe fields in an atomic system and induce a charge ex- 



citation of the write port and a phase excitation of the 
read port, respectively. This hamiltonian has been writ- 
ten supposing that the offset gate charge and loop flux 
have been compensated to operate at the optimal bias 
point where the charge ^ and the flux ^ have zero 
mean value in both the ground |0) and first excited |1} 
states of Hbox- Under these conditions, the qubit is min- 
imally sensitive to charge and flux noise 0. 

If we keep these two lowest states in the Hilbert space 
of the circuit we obtain an effective qubit hamiltonian 
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The photon annihilation operator a is related to S by 

J a + 



{Ef/2Eg) 



1/4 



which represents the decomposition of the gauge- 
invariant phase difference into annihilation and creation 
operators of the "plasma" mode whose bare frequency is 
LUp. The operators ax and az are the Pauli spin opera- 
tors and E^ is the single electron charging energy of the 
readout junction. In this effective hamiltonian, the ex- 
pansion of cos 6 is carried out only to the flrst anharmonic 
term, which describes the non-linear resonator dynamics 
with sufficient accuracy for a bifurcation readout. 

Let us describe the role of each term in (1). The first 
term describes the infiuence on the qubit of the charge 
port drive which is used to manipulate its state. The 
second term is the Larmor term ujqi — Ej/h. We have 
supposed here that the ratio Ej /Ec is sufficiently small 
that corrections to the Larmor frequency involving Ec 
are small. To model the behavior of qubit samples with 
an appreciable Ej/Ec ratio, we would keep higher order 
terms, yielding renormalized values of the coefficients in 
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(1). The third term describes the dominant couphng be- 
tween the qubit and the resonator. Note that this term 
commutes with the hamiltonian of the qubit when U = 0, 
offering the possibihty of quantum non-demohtion mea- 
surements. The fourth term describes a decrease in the 
frequency of the resonator when its photon population 
increases Finally, the fifth term describes the ex- 

citation of the resonator by the drive current applied 
through the phase port. When the drive current is in- 
creased while its frequency is sufficiently below cUp the 
system becomes metastable with two possible dynamical 
states with different oscillation amplitudes, i.e. two pos- 
sible photon populations . We exploit this bistability 
for our readout, which we describe in the next section. 

III. QUBIT READOUT 

It is clear from the hamiltonian (1) above that the dy- 
namics of the non-linear resonator depend on the value 
CTz = ±1 corresponding to the state of the qubit. In par- 
ticular, the small oscillation "plasma" frequency oj^^ = 
u)p (1 ± A) varies with the qubit state. We probe the non- 
linear resonator by sending down the phase port trans- 
mission line a microwave pulse with carrier frequency 
Lu — ujp ~ Alu, such that the detuning Auj > ^(^p 
where Q is the quality factor of the plasma resonance 
|l7j |. In our circuit, the damping of the plasma res- 
onance arises from the characteristic transmission line 
impedance Zc = and thus Q = ZcCujp ~ 10 — 20. 
For this value of detuning, when ramping up the drive 
current the resonator switches from one dynamical 
state to another when 



smaller than Is [uj,ujp{l — A)], is still higher than the re- 
verse bifurcation "retrapping" current I§. This latching 
property conserves the information about the qubit state 
acquired during a small time interval in the resonator 
and allows us to probe the reflected phase (p during a 
time typically longer than r^. 

In Fig. 2, we present a typical histogram of the re- 
flected drive signal phase (f) corresponding to a drive cur- 
rent Irf which causes the resonator to switch, on average, 
half of the time. The histogram has 800,000 counts ac- 
quired in 200 ms. For qubit measurements shown later, 
histograms with only 10,000 are used. The shape of the 
readout pulse used is schematically shown in the inset of 
Fig. 2. The rise time of the pulse is typically 20 — 40 ns, 
the maximum current I^j'^^ is applied for 40 — 120 ns 
and the latched section lasts 120 ns, during which the 
recorded reflected signal phase (j) is bimodal, with values 
differing by about 124°. We have chosen the phase refer- 
ence so that the value — ^low = —62° corresponds to 
the resonator in its initial state, while (j) — 0high = 62° 
corresponds to the resonator having switched. We define 

the switching probabihty P^witch {^I^f"", Aw, \az\ , 
where j^*} is the state of the qubit, as the weight of the 
histogram that lies above (f> — — Q. 

IV. COHERENCE RESULTS 

We now present experimental results on two different 
qubit samples whose characteristic parameters are listed 
in Table I, along with a summary of our results. In the 
figures that follow, we only show data for sample A. 



Irf>lB{c0,L0f) 

where Ib is the bifurcation current with expressions given 
in T^. Therefore, by choosing the maximum pulse am- 
plitude 



Ib [lo, LOp (1 - A)] < I^^^ < Ib [lo, cop (1 + A)] 

we can determine, by measuring if the resonator has 
switched or not, whether the qubit was in state |0) or 

The dynamical states of the resonator differ in both 
the amplitude and phase of the forced oscillations at fre- 
quency u. In this work, we have chosen to use a reflec- 
tometry setup in which all the information about the res- 
onator state is carried by the reflected drive signal phase 
(j). This last property occurs because the probed circuit is 
not intrinsically dissipative (in absence of quasi-particles, 
which is very well realized in our measurements) and 
the power reflected from the chip is equal to the inci- 
dent power in steady state. A further advantage of our 
non-linear resonator is that the switching is strongly hys- 
teretic. Once a switching event has occurred we can de- 
crease the drive current Irf to a value which, while much 
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FIG. 2: Typical histogram of the phase of the reflected signal 
in the JBA readout when the maximum rf drive current is 
chosen so that the resonator switches approximately half of 
the time. The switching probability Pswitch is deflned as the 
fraction of the histogram lying above = 0. The inset shows 
schematically the envelope of the readout pulse sent to the 
phase port. The qubit influences the switching probability 
during the time interval Tm which here was 40 ns. 
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TABLE L Parameters for two measured qubit samples. The readout frequency was 1.55 GHz and 1.70 GHz for samples A and 
B, respectively. The detuning was 6% of tOp. The parameter rj is the discrimination power of the readout. 



Sample 


tJoi/27r(GHz) 


Ej/Ec 


21, typical (/is) 


T2(ns) 


recho (ns) 






^cxp / ^calc 


A 


9.513 


2.7 


4.0 


320 


400 - 500 


.48 


.70 ± .05 


.69 


B 


18.989 


6.0 


1.0 


300 


300 


.61 


.70 ± .05 


.87 



We first characterized our readout by measuring 
^switch as a function I^p'^ and as shown in Fig. 3. 
The blue circles correspond to the qubit in its ground 
state, obtained by letting the qubit relax spontaneously, 
while the red circles correspond to the qubit in its first 
excited state obtained by applying a tt pulse, which will 
be discussed below. An important remark is that only 

a slight change in shape of Pswitch ij^f'^^ between the 
two qubit states is observed, which indicates that the 
switching process itself does not contribute strongly to 
the relaxation of the qubit. In cases where the readout is 
suspected to induce significant relaxation, the switching 
probability curve for the qubit excited state displays a 
pronounced kink and can be obtained by a weighted av- 
erage of the observed curve for the ground state and the 
prediction for the excited state 0, llOl • The discrimina- 
tion power of the qubit readout is defined as 
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FIG. 3: Switching probability as a function of maximum drive 
current and qubit state for sample A. Dotted line represents 
value of drive current at which maximal discrimination power 
is observed. The width in current of the curves is in reasonable 
agreement with numerical simulations (data not shown) . The 
solid line connects the observed data points in the |0) state 
and the dashed line is a copy of the solid line horizontally 
shifted to overlap the jl) state data at low values of Pswuch- 



T] = max [Pswitch {{o■z)^s, = 1) - Pswitch {{<^z)^s, = -1)] 

and its observed {rjexp) and predicted (77 raic) values are 
given in Table I. Numerical simulations |2l| of the full 
circuit have been used to compute the predicted values 
of 77. Note that several competing factors enter this cal- 
culation, yielding similar values for samples A and B. The 
error bars reflect uncertainties in the values of stray re- 
actances on chip and the precise resonator temperature. 

The observed discrimination power is about 15 — 30 % 
smaller than expected, and we attribute this loss to spu- 
rious on-chip defects. In a set of experiments to be de- 
scribed in a later publication, we used two readout pulses 
in succession to determine that a 15 — 30 % loss of qubit 
population occurs even before the resonator is energized 
to its operating point. As photons are injected into the 
resonator, the effective qubit frequency is lowered due to 
a Stark shift via the phase port [23|- When the Stark 
shifted frequency coincides with the frequency of an on- 
chip defect, a relaxation of the qubit occurs. Typically, 
the qubit frequency spans 200 — 300 MHz before the state 
of the qubit is registered by the readout, and 3 — 4 spu- 
rious resonances are encountered in this range. 

For future measurements, we have developed a new 
method to counter this effect. When applying a read- 
out pulse via the phase port, we apply a compensating 
pulse via the charge port which Stark shifts the qubit 
to higher frequencies. When balancing these pulses, we 
have successfully reduced the net frequency shift to below 
20 MHz and have minimized population loss to defects 
before the resonator switches. To increase the expected 
discrimination power to unity, we must use samples with 
either a larger qubit Ej or a stronger phase coupling be- 
tween the qubit and readout resonator. The latter can 
be accomplished by using a resonator with two Josephson 
junctions in series. 

Having characterized our readout discrimination 
power, we performed a series of experiments to assess 
the coherence of our qubit, namely the measurements of 
Ti, T2, Tccho and T2. These times characterize the decay 
of the excited state population after a tt pulse, the de- 
cay of Ramsey fringes, the decay of the echo signal after 
a (7r,7r/2, tt) pulse sequence and the decay of the Rabi 
oscillations, respectively. 

We first applied to the charge port a pulse at the Lar- 
mor frequency woi of varying duration r and amplitude 
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FIG. 4: a) Rabi oscillations of the switching probability 
of qubit sample A as a function of the duration t of a 
square pulse applied on the gate with maximum amplitude 
U^^" = 0.12mV. Solid green curve is an exponentially decay- 
ing sinusoidal fit with T2 — 1.6 us. Total acquisition time is 
3mn and duty cycle is 16/is, set by Ti (see below), b) Rabi 
oscillation frequency measured in (a) as a function of C/^^^". 
Green line is expected linear dependence. 



U"j:'^^, which performs a rotation of the qubit about ax, 
followed by a readout pulse on the phase port. The result- 
ing Rabi oscillations in the switching probability signal 
are plotted in Fig. 4a for varying r and fixed U^^^. Near 
T = we observe the Pswitch corresponding to qubit be- 
ing in the |0) state. As the pulse length increases, Pswitch 
increases, goes through a maximum where the qubit is 
purely in the |1) state, defining at this point the length 
of a TT pulse. The switching probability then decreases 
back to the |0) state value, indicating a full 2n rotation 
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FIG. 5: Decay of the excited state switching probability after 
preparing the qubit in the excited state by a tt pulse, as a 
function of the waiting time between the preparation pulse 
and the readout pulse. Data for sample A. Solid green curve 
is an exponential fit with a 3.2/xs decay constant. 



of the Bloch vector. This pattern repeats itself but with 
diminishing contrast. The decay time T2 is in the range 
0.8 — 1.7 /j-s depending on the sample and precise biasing 
condition. The Rabi oscillation frequency lor is plotted 
as a function of U^^^ in Fig. 4b. A linear dependence 
of u}R with U^j^'^ is observed, in agreement with theory. 
The shortest tt pulse we generated was 2 ns long, and was 
used in the echo experiments described below. 

Having calibrated the tt pulse, we then performed a 
qubit energy relaxation measurement by introducing a 
waiting time between the tt pulse and the readout 
pulse. The decay of Pswitch with t^j, shown in Fig. 5, 
is well fitted by a single exponential, defining Ti. For 
sample A, Ti was in the range 2.5 — 5 /iS, and for sample 
B, Ti was between 1.0 — 1.3 /iS. The values of Ti ob- 
tained with our dispersive readout are comparable with 
the results of Vion et al. , and are significantly shorter 
than the values expected from coupling to a well ther- 
malized 50 VL microwave environment shunting the qubit. 
The loss mechanisms giving rise to the observed energy 
relaxation are not understood at this time. 



Following measurements of the qubit energy relax- 
ation, we performed a Ramsey fringe experiment to de- 
termine the phase coherence of the qubit. In this exper- 
iment, two 7r/2 pulses were applied to the charge port 
of the qubit at a frequency 10 — 20 MHz detuned from 
ciJoi followed by a readout pulse on the phase port. A 
free evolution time At was introduced between the two 
7r/2 pulses. In Fig. 6, Pswitch is plotted as a function of 
/S.t. In the Ramsey sequence, the first 7r/2 pulse tips the 
Bloch vector from the north pole to the equatorial plane. 
During the time At, the Bloch vector precesses around 
the equatorial plane and is then rotated again by the sec- 
ond 7r/2 pulse. For At = 0, the two 7r/2 pulses back to 
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FIG. 6: Ramsey fringes obtained with two 7r/2 pulses sep- 
arated by the time interval At. The pulse frequency was 
detuned from the Larmor frequency by 20MHz. The green 
curve is a exponentially decaying sinusoid fit. The decay time 
T2 is 320ns. Same acquisition conditions as in Fig. 4. 

back act as a single tt pulse and the observed value of 
-Pswitch corresponds to the qubit being in the |1) state. 
As At increases, ^switch decreases until it reaches the 
value corresponding to the qubit being in the |0) state, 
corresponding to a free evolution time At in which the 
Bloch vector makes a tt rotation in the equatorial plane. 
The switching probability then continues to increase for 
larger values of At until it reaches a maximum value, cor- 
responding to a time At where the Bloch vector makes 
a full 2tt rotation in the equatorial plane. This oscilla- 
tory pattern then repeats but with decreasing contrast 
corresponding to the loss of phase coherence with time. 
The Ramsey fringes decay in a time T2 which has a com- 
ponent due to energy relation and one due to pure de- 
phasing: I/T2 = l/(2Ti) + l/T^, where represents 
pure dephasing. In our measurements, T2 is dominated 
by pure dephasing. For sample A, T2 = 320 ns, and for 
sample B, T2 = 300 ns. 

In order to correct dephasing of the qubit due to low 
frequency noise [^.E^, we performed an echo experiment 
in which we inserted a tt pulse in the middle of the two 
7r/2 pulses of the Ramsey sequence. A set of Ramsey 
fringes and its corresponding echo decay are shown in Fig. 
7 for sample A. For this sample, the decay constant was 
increased to 400 — 500 ns using the echo technique. For 
sample B, the echo technique did not increase the phase 
coherence time. We believe that for sample B, which has 
a large ratio of Ej/Ec and is protected from 1// offset 
charge noise, the dominant source of dephasing is due 
to broadband noise emanating from residual photons in 
our readout resonator , thus explaining the inefficacy 
of the echo sequence. It is possible that the 50 envi- 
ronment shunting the qubit on the phase port side was 




At (ns) 

FIG. 7: Result of an echo experiment where a tt pulse was kept 
in the middle of the two n/2 pulses separated by interval At 
(black dots) . The corresponding Ramsey fringe data is shown 
with red circles. The thick black curve is an exponentially 
decaying fit. 

not fully thermalized to the refrigerator temperature of 
10 mK. For sample A, where an improvement was ob- 
served with the echo sequence, there are likely two con- 
tributing factors. First, the ratio Ej/Ec is much smaller 
and offset charge noise played a stronger role. The low 
frequency component of this noise can be corrected us- 
ing an echo sequence. Second, we added more cryogenic 
attenuation in the transmission lines directly coupling 
to the phase port to reduce the resonator temperature, 
thereby potentially reducing the number of excess pho- 
tons in the readout resonator and their associated de- 
phasing. 



V. CONCLUSION 

In conclusion, we have successfully implemented a non- 
linear dispersive readout of the quantronium qubit using 
the Josephson Bifurcation Amplifier. The readout speed 
and discrimination power show a significant improvement 
when compared with the DC switching readout used in 
the original quantronium measurements Q. Perhaps 
even more important, in the present readout scheme, the 
total measurement time is much smaller than Ti , and it is 
possible to carry out experiments with multiple readout 
pulses to determine the information flow during a qubit 
readout and to account for any losses in qubit popula- 
tion. This important aspect can be used to determine 
the degree to which the measurement is quantum non- 
demolishing, and will be described in later publications. 
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